The 17/5 spectrum of the Kelvin-wave cascade 
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Direct numeric simulation of the Biot-Savart equation readily resolves the 17/5 spectrum of the 
Kelvin-wave cascade from the 11/3 spectrum of the non-local (in the wavenumber space) cascade 
scenario by L'vov and Nazarenko. This result is a clear-cut visualisation of the unphysical nature 
of the 11/3 solution, which was established earlier on the grounds of symmetry. 
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A cascade of Kelvin waves (kelvons) — distortion waves 
on quantized vortex filaments — plays a crucial role in the 
decay of superfluid turbulence at T = by providing 
the mechanism of energy transfer [ll-[3 to shorter length 
scales, where the dissipation due to phonon emission be- 
comes efficient [5|, Q (for a review of the theory of su- 
perfluid turbulence in the T — )► limit, see Ref. Q). 
Our previous analytic 0) and numeric Q studies of the 
Kelvin-wave cascade in the regime of weak turbulence 
revealed the spectrum 
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for the kelvon occupation numbers. Since then this 
result — implying the locality of the three-kelvon elas- 
tic collisions in the wave number space — was consid- 
ered well-established. However, recently, Laurie, L'vov, 
Nazarenko, and Rudenko claimed the three-kelvon (a.k.a. 
six- wave) collisions to be essentially non-local 8]. An al- 
ternative theory of the Kelvin- wave cascade was then put 
forward by L'vov and Nazarenko [9[, with the spectrum 
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Our subsequent analysis of the problem [10] brought us 
to the conclusion that the locality of the Kelvin- wave cas- 
cade can be rigorously proven on the basis of geometric 
symmetries alone. Moreover, we argued that the spec- 
trum ^ is characteristic of (uncontrollable) approxima- 
tions, as well as errors, violating the so-called tilt sym- 
metry. We have also remarked that proposed in Ref. [8!| 
'additional' terms of the collisional integral — absent in 
our analysis [3| — are likely to be the source of the mis- 
take. The reaction of the community to the argument 
of Ref. [iO] remains essentially reserved [ll[. Therefore, 
we find it instructive to visualize the argument by first- 
principle direct simulations. 

Vortex dynamics at T = is described by the (asymp- 
totically exact in the hydrodynamic limit) Biot-Savart 
equation (BSE). The BSE is notoriously difficult to inte- 
grate numerically due to the essentially non-local (in real 
space) interaction between vortex-line elements. More- 
over, it is this interaction that accounts for kelvon kinet- 
ics in the weak-turbulence regime in question. To effi- 
ciently simulate the BSE, in our Ref. [Jl, we developed 



the so-called scale-separation scheme, which reduces the 
simulation cost of BSE to that of a local model. Back 
in 2005, the efficiency of the scheme was sufficient to 
unambiguously resolve the difference between the spec- 
trum k~^'^'^ and the previously conjectured k~^, using 
a standard desktop computer. Given that the difference 
between 11/3 and 17/5 is of the same order as the differ- 
ence between 17/5 and 3, it is expected that the numeric 
protocol of Ref. f4| should allow one to readily resolve 
17/5 from 11/3 with a modern laptop. 




FIG. 1. (Color online) The time evolution of the Kelvin-wave 
spectrum plotted as Ukk^'' ■ The initial condition (f = 0) 
corresponds to the spectrum ([2]). As the time t (measured 
in the units of the period of the fastest mode in the system) 
progresses, the spectrum ([2]) is transformed by a wave prop- 
agating from high wavenumbers to the smaller ones into the 
expected n^ oc fc~^^'^ spectrum. 

The results of our simulation are presented in Fig. [T] 
For convenience of visual inspection we plot n^fc^^/^ 
so that the spectrum ([1]) corresponds to a horizontal 
line. We use the same routine as in Ref. 4], with the 
only difference that the initial condition (labelled with 
t — Q) now corresponds to the spectrum ^ rather than 
rife ex k~^ . We show the evolution of the spectrum at two 
instants, i = 4 and t — 600 (in the units of the period of 
the fastest mode in the system). The spectra follow from 



time averaging of the wavenumber distribution over in- 
tervals At <t suggested by the quasi-steady-state nature 
of a cascade. In order to filter the inherent noise between 
neighbouring wavenumbers, we obtain rik after averaging 
over the wavenumber range [1.5~^"fc, 1.5^'^fc]. The sim- 
ulation took 12 hours on a laptop with a 2.4 GHz proces- 
sor. We clearly see that the distribution nt oc fc~^^/^ is 
being re-structured into Uk oc fc~^^'^, the transient hav- 
ing the form of a wave propagating from large to smaller 
wavenumbers, in a direct analogy with the picture seen in 
Ref . [J] . It is important to note that during the transient 
the amplitude of the fc~^^/^ tail (the height of the hori- 
zontal piece of the curve) increases. This is a qualitative 
feature immediately ruling out the spectrum ([2|). 



In conclusion, we make two remarks of a more general 
character. First, we see that there is no need for de- 
veloping approximate local (in real space) models since, 
in numerical simulations of kelvons on individual vortex 
lines, the full spatially non-local model can be efficiently 
solved at the expense of a local one (and when the cou- 
pling between separate lines is important, no local model 
is applicable in principle). Second, the efficiency of the 
scale-separation scheme calls for the implementation of 
tree codes in simulations of vortex tangles. 
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